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A precision method for determining 
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V u b\ using the full range in q of B 
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iriu data is presented. 

At large q A the form factor is taken from unquenched lattice QCD, at q 2 — we impose a model 
independent constraint obtained from B — > tttt using the soft-collinear effective theory, and the shape 
is constrained using QCD dispersion relations. We find \V ub \ = (3.54 ± 0.17 ± 0.44) x 10" 3 . With 
5% experimental error and 12% theory error, this is competitive with inclusive methods. Theory 
error is dominated by the input points, with negligible uncertainty from the dispersion relations. 



The remarkable success of the B-factories have lead to 
a new era for precision results in the CKM sector of the 
standard model. For \V u b\ inclusive and exclusive mea- 
surements from semileptonic decays should yield a precise 
value, but must surmount the now dominant theoreti- 
cal uncertainties. For inclusive decays measuring \ V u b\ is 
more difficult than \V c b\ because cuts make observables 
either sensitive to a structure function which demands 
input from radiative decays, or require neutrino recon- 
struction. The heavy flavor averaging group (HFAG)'s 
average from inclusive decays based on operator product 
expansion techniques is 10 3 |T4b| = 4.7±0.4 1]. Exclusive 
techniques for |U c f,| use heavy quark symmetry (HQS) to 
normalize the form factors. For \ V u b\ from B — > t:£9 sym- 
metry techniques fall flat, and model independent form 
factor information relies on precision lattice QCD. 

Recently, the Fermilab |2| and HPQCD 3] groups have 
presented unquenched lattice results for B — ► tt form 
factors. Uncertainties in the discretization restrict the 
kinematics to pions that are not too energetic E n < 
1 GeV, which for the invariant mass of the lepton pair 
is 15 GeV 2 < q 2 < 26.4 GeV 2 . Unfortunately, since the 
phase space goes as |p w | 3 , there are less events and more 



experimental uncertainty in this region. For B 
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dT/dq 2 = N \V ub \ 2 \f+(q 2 )\ 
where N = G 2 F / (24w 3 ) . For example, Belle Q found 
" 3.87± 0.70 ± 0.22±g-|f (FNAL) 



(1) 



10 |Kb|,j2> 16 



4.73 ± 0.85 ± 0.27^50 (HPQCD) 



(2) 



where the errors are statistical, systematic, and theoret- 
ical. In quadrature this is an uncertainty of ~ 25%. 
The latest Babar, CLEO, and Belle average is ||, 



Br(S° 



r + t~v) = (1.39 ±0.12) x 10" 



(3) 



which should yield \ V u b\ at the ~ 5% level. So far extrac- 
tions of \V u b\ from the total Br rely on QCD sum rules 6] 
and quark models for input. For example, HFAG reports 
results on Bi(B — > {it, p, io\lv) that lead to central values 



10 3 |K b | = 2.9 to 3.9 1]. Due to the uncertainty they do 
not currently average over exclusive extractions of \V u b\- 

In this letter we present a model independent exclu- 
sive method for determining the entire B — > 7r form fac- 
tor f+(q 2 ) and thus \V u b\- A total uncertainty 6\V u b\ — 
13% is achieved by combining 1) the unquenched lat- 
tice results 0, j 2) a constraint at q 2 = derived 
from SCET jjj and B — > 7T7t data, which determines 
|Kif,|/+(0), and 3) dispersion relations and analyticity 
which allow us to interpolate over the entire region 
of q 2 by bounding the shape of f+(q 2 ) between input 
points [alii- The SCET constraint induces an additional 
implicit functional dependence on | V^t. | in the form fac- 
tors. Our first analysis uses just the total Br, yielding an 
analytic formula for \ V u b\. The second includes q 2 -spectra 
with a x 2 minimization which allows the experimental 
data to constrain the theoretical uncertainty. A differ- 
ent approach for including the g 2 -spectra was developed 
in |l£| based on the Lellouch distribution method 

Analyticity Bounds. We briefly review how analytic- 
ity constrains the B — ► 7r form factors, /+ and /o, re- 
ferringto 0, 0, for more detail. Our notation fol- 
lows [13], and we set t± = (m_e ± m-*) 2 ■ Suitable mo- 
ments of a time ordered product of currents, n^g 2 ) = 
i Jd i xe lqx (Q\TJ' l (x),P u {Q)\Q) can be computed with an 
OPE in QCD and are related by a dispersion relation to 
a positive definite sum over exclusive states 



Im II^= /[p.s.] 6 (q- P B*) (0 1 J ]v I Bit) (Btt\ J» \ 0) 
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Keeping this first term bounds a weighted integral over 
t + < t < 00 of the squared Bit production form factor. 
Using analyticity and crossing symmetry this constrains 
the shape in t = q 2 of the form factors for B — > tt in the 
physical region < t < t_. The results are simple to 
express by writing each of f+(t), fo{t) as a series 



/(*) 



1 



P{t)4>(t,t ) 



J2a k (to)z(t,t ) k , (5) 



fc=0 



with coefficients at that parameterize different allowed 
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functional forms. The variable 



z(t,to) 



t-Vu 



to 



VU=t + - t 



(6) 



maps t + < t < oo onto \ z\ = 1 and — oo < t < t + onto z E 
[—1,1]. to is a free parameter that can be chosen to attain 
the tightest possible bounds, and it defines z(t 0l to) = 0. 
We take t = 0.65 1- giving -0.34 < z < 0.22 for the 
B — > 7r range. In Eq. © the "Blaschke" factor P(t) 
eliminates sub-threshold poles, so P(t) = 1 for /o, while 
P(t) = z(t;m 2 B «) for / + due to the B* pole. Finally, the 
"outer" function is given by 
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(*+-*) 



(a+l)/4 



(t+^to) 1 / 4 



-(6+3) 



(y/U^i+y/U-t-) 



a/2 



(7) 



where n.j = 3/2 and for / + : (A" = 487r, a = 3, b = 2), 
while for f : (K = 16n/(t+t-),a = 1,6 = 1). Here xT 
corresponds to the lowest moment of n(q 2 ) computed 
with an OPE. At two loops in terms of the pole mass 
and condensates and taking /i = m& 0, El 



(o) _ 3 [1 + 1.140 a s (mb)] m b (uu) 
x f+ _ " 



32n 2 ml 



(a s G 2 ) 



*f< 



(0) 



[1+0.751 a s (m&)] m b (uu) (a s G 2 ) 



with mf,(uu) 

pole 



8tt 2 



-0.076 GeV 4 , (a s G 2 ) 



127rm 4 ' 



(8) 



0.063GeV 4 . We 
4.88 GeV as a central value. With Eq. © 
the dispersive bound gives a constraint on the coefficients 



riA 

£« 2 fc <i 

fe=0 



(9) 



for any choice of ha- 

Eqs. (J5J and 10 give only a weak constraint on the 
normalization of the form factor /+. In particular, data 
favors ao ~ 0.02, so a 2 , <C 1. The main power of ana- 
lyticity is that if we fix f + (q 2 ) at ua input points then 
it constrains the q 2 shape between these points. With 
Ua = 5 the error from the bounds is negligibly small rel- 
ative to other uncertainties, as we see below (our analysis 
is also insensitive to the exact values of Xj or m b)- The 
bounds can be strengthened using heavy quark symme- 
try or higher moments of n(<? 2 ) |l2j, but since this uncer- 
tainty is very small we do not use these improvements. 

Input Points. A constraint at q 2 — is useful in pin- 
ning down the form factor in the small q 2 region. Here we 
implement a constraint at q 2 = on |I4b|/+(0) that fol- 
lows from a B — ► tttt factorization theorem derived with 
SCET @ ■ The result holds in QCD and uses isospin sym- 
metry and data to eliminate effects due to the relative 
magnitude and strong phase of penguin contributions. 



Manipulating formulas in we can write the result in 
terms of observables 



|Kf,|/+(0) = 



64tt Br(B~ -> tt q tt-) 



ml ft r B -\V ud \ 2 G 2 F 
(Ci + C 2 )t c - c 2 



1/2 



(10) 



cl-c 2 



1 + (a s (m b ), 



A 



QCD 



where C\ — 1.08 and C 2 = —0.177 are parameters in the 
electroweak Hamiltonian at /i = m b (we drop the tiny 
and t c is a hadronic parameter whose deviation 
from 1 measures the size of color suppressed amplitudes. 
In terms of the angles (3, 7 of the unitarity triangle and 
CP-asymmetries S v + V - and C„+„- in B — ► 7r + 7r~, 



u (i+^+jr- cos2j3 + S v+v - sin 2/3) 
V 2 sin 7 

withR c =[Br(B° -> ir+it-)T B -]/[2Br(B- -> 7r 7r - )TB°], 
and JW = (1- C 2 + ^_ -5 2 +w _)V2. Eqs . im _ 

prove on relations between B — > 7T7r and _B — > 7r^^ 
derived earlier, such as in Ref. 0], because they do not 
rely on expanding in a s (\/m A) or require the use of 
QCD sum rules for input parameters to calculate t c . 
Using the latest B — ► hit data yj, Eq. I|10|l gives 



/f n = |Kfe|/+(0) = (7.2±1.8)xlO- 



(12) 



This estimate of 25% uncertainty accounts for the 10% 
experimental uncertainty, and ~ 20% theory uncertainty 
from perturbative and power corrections. The experi- 
mental uncertainty includes 7 = 70° ± 15° which cov- 
ers the range from global fits and that preferred by the 
SCET based B — > irn method from Ref. [lj|. As noted 
in the dependence of |V^(,|/+(0) on 7 is mild for larger 
7's. Estimates for perturbative and power corrections to 
Eq. (|10|) are each at the ~ 10% level even when "chirally 
enhanced" terms are included [14Lll6| . 

Next we consider lattice QCD input points, which 
are crucial in fixing the form factor normalization. Tech- 
nically, using staggered fermions might add model depen- 
dence from the (detM) 1 / 4 trick. We take the remark- 
able agreement in as an indication that this model 
dependence is small. Using the unquenched MILC config- 
urations, Refs 0j0 find consistent results with different 
heavy quark actions. As our default we use the Fermilab 
results since they have a point at larger q 2 : 



fl n = /+(15.87) = 0.799 ± 0.058 ± 0.088 , 
f 2 n = /+ (18.58) = 1.128 ± 0.086 ± 0.124 , 
/f n = /+ (24.09) = 3.262 ± 0.324 ± 0.359 . 



(13) 



The first errors in (|13|l are statistical, ±cr,-, and the second 
are 11% systematic errors, ±yf? n , with y = 0.11. For the 
lattice error matrix, we use Eij = of J,j + y 2 f\ n f( n , which 
takes Oi uncorrelated and includes 100% correlation in 



3 



the systematic error. Of the eleven reported lattice points 
we use only three at separated q 2 . This maximizes the 
shape information while minimizing additional correla- 
tions that may occur in neighboring points, for example 
from the chiral extrapolation. 

Chiral perturbation theory (ChPT) gives model inde- 
pendent input for / + (and /o) when E„ ~ m n , namely 



f + {q 2 (E 7T )) = 



gfBmB 



1- 



-(f) 



(14) 



where g is the B*Btt coupling and fs the decay con- 
stant. Possible pole contributions from the low lying 
J n = + ,l + ,2 + states vanish by parity and angular 
momentum conservation. The first corrections scale as 
E^/A, where A ~ 600 MeV is the mass splitting to the 
first radially excited 1~ state above the B*. We take g = 
0.5. This is compatible with D* decays using heavy quark 
symmetry. Updating the ChPT fit in by including 
both T(D* + ) and D* Br-ratios, gives gn*D^ — 0.51 (at 
an order where there are no counterterm operators and 
no l/m c corrections absorbed in g). For the lattice aver- 
age Hashimoto gives f B = 189 MeV. Thus, 



/£ = /+ (26.42) = 10.38 ±3.63, 



(15) 



where this fairly conservative 35% error is from uncer- 
tainty in gfs, and from the m T /A ~ 23% corrections. 
Determining f + . To determine f+(t) we drop afc>6 in 



) 1 ' 2 which properly 
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Eq.JoJ, and take a*, — > a*,{l - 

bounds the truncation error 20]. The / 0_ * input points 
then fix ao-4 as functions of 05. Functions that bound 
f+(t) are determined from the maximum and minimum 
values of 05 satisfying © with tia = 5. Thus we solve 

18.3ao±3.96ai±0.857a 2 ±0.185a3±0.0401a 4 (16) 

+0.00887a 5 = f/\V ub \ , 
37.8ao±0.960ai±0.0244a 2 ±0.000619a 3 ±1.57xl0~ 5 a4 

±4.00xl0~ 7 a 5 = f 1 , .. . , 
304.0oo-103.6oi+35.3a 2 -12.0o3+4.10o4-1.49a 5 = / 4 , 

a + a l + a 2 + a 3 + a 4 + a 5 = 1 • 

In Eq.JSJ this yields two solutions, F±, with parameters 

f+(t) =F ± (t,{h/\V ub \,h,h,h,U}). (17) 

To see how well these solutions bound the form factor 
we fix \V u b\ = 3.6 x 10 -3 , /* = ff n and plot the bounds 
as the two black solid lines in Fig. The curves lie on 
top of each other. For comparison we show dashed lines 
for the bounds on / + and /o obtained using four lattice 
points (shown as dots). With these inputs the constraint 
/+(0) = /o(0) is less effective than using the SCET point. 

\V U b\ from total Br -fraction. Equating Eq.© with the 
theoretical rate obtained using Eas.p7jl gives an analytic 
equation for \V u b\. With f % = f? n the solution is 
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FIG. 1: Upper and lower bounds on the form factors from 
dispersion relations, where q 2 = q 2 /m%* and the (1— q 2 ) factor 
removes the B* pole. The overlapping solid black lines are 
bounds F± derived with the SCET point, 3 lattice points, 
and the ChPT point (diamonds with error bars). The dashed 
lines are the bounds derived using instead four lattice points 
(shown by the dots). Input point errors are not included in 
these lines, and are analyzed in the text. 



Type of Error 


Variation From 


<5|I4i>| r 


5\V ub f 


Input Points 
Bounds 

pole 

ml 
OPE order 


l-cr correlated errors 
F+ versus F- 

4.88 ± 0.40 
2 loop — > 1 loop 


±14% 
±0.6% 
±0.1% 
-0.2% 


±12% 
±0.04% 
±0.2% 
±0.3% 



(4.13 ±0.21 ±0.58) x 10" 



(18) 



TABLE I: Summary of theoretical uncertainties on | V u b I • Re- 
sults are shown for an analysis from the total branching frac- 
tion, <5|Vui,| Br , and from using the dF/dq 2 spectrum, ^IV^I 9 . 
For the input point error we quote the average from F±. 

The first error is experimental, 5.2%, propagated from 
Eq.©. The second error, 14%, is from theory and is 
broken down in Table [fl It is dominated by the input 
points. The bound uncertainty from the choice of so- 
lution is < 1% (but would grow to ±12% without the 
SCET point). The error from mb and the order in the 
OPE and are very small because shifts in the normaliza- 
tion through xfl are compensated by shifts in the a n co- 
efficients, except for the last term 05 which gives a small 
contribution. To ensure consistency with the dispersion 
bounds the input point uncertainty is calculated using 
the Lellouch-method of generating random points from 
Gaussians 11|, giving 10 3 |Kb| = (3.96 ± 0.2 ± 0.56) . 
Our distributions were determined using Eqs. (|12I13I15|) 
and the correlation matrix E^. Taken individually the 
SCET and ChPT points give ~ 5% error, so the lattice 
uncertainty dominates. 

I Vub I from q 2 spectra. Results for partial branching 
fractions, (Br° xp ± <5BrU, over different bins in q 2 are 
also available. Cleo jij and Belle present results 
for 3 bins with untagged and ir + semileptonic tags re- 
spectively. Babar recently presented total rates from 
hadronic & leptonic ir + and ir° tags as well as ir + semilep- 
tonic tagged data in 3-bins and untagged data over 5- 
bins. By fitting to these 17 pieces of data with Minuit 
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FIG. 2: Results from the \ 2 fit of \V ub \ and f°~ 4 to the q 
spectra (q 2 — q 2 /m%»). The two solid lines are obtained using 
either the F+ or F- solutions from Eq. 117H . The two dashed 
lines repeat this analysis without using the SCET point. 
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FIG. 3: The curves are as in Fig.2, but for the decay rate. 



we exploit the q 2 shape information. To do this we define 



E 



[B^-B Ti (V ub ,F ± )f , [/°-/°] 



012 



(Sf 



412 3 



(Sf ) 2 



(19) 



4^2 



E [ftn-fM-fUE-'h 



and minimize \ 2 as a function of \V ub \ and J 0-4 , x 2 
contains both experimental and theoretical errors, with 
-E -1 the inverse error matrix. By allowing /°~ 4 in F± to 
move away from f?~ A the theoretical rate is allowed to 
adjust itself based on the experimental q 2 shape. 
Minimizing i|19|) gives X 2 /(dof) = 1.04 and 



V„ h \ = 



(3.54 ±0.47) 



x 10 



-3 



(20) 



Results for f+(q 2 ) and dT/dq 2 are shown by the black 
solid curves in Figs. and Eq. (J2DJ) has a total error of 
13%. If we fix /°~ 4 = /°~ 4 then the experimental error is 
4.9%, ie. S\V ub \ = ±0.17. The remainder, S\V ub \ = ±0.44 
is from the input points, so the q 2 spectra brought this 
theory error down to 12%. Other uncertainties are small 
as shown in Table [I] The experimental spectra favor a 
larger form factor between the lattice and SCET points. 
This decreases the value of \V ub \ from that in l|18|) . Using 
Eqs. (tT2TPol) this fit yields 

/+(0) = 0.227 ±0.047, g f B = 96 ± 29 MeV , (21) 

consistent with our inputs. This /+(0) has 21% error. 



If we entirely remove the SCET point f° from Ea. i|19|) 
then we obtain a fit that uses only semileptonic data, 
shown by the dashed red lines in Figs. and The 
spectrum is now determined less precisely at small q 2 , 
since this data only bounds the area in the smallest q 2 - 
bin. The result is \V ub \ = (3.56 ± 0.48) x 10~ 3 . It has 
the same input point error as Eq.J3U|) and a somewhat 
larger bound error, S\V ub \ = 1.8%. Turning the use of 
Ea. (fT2*| around, we can combine it with /+(0) to get an 
independent method of fixing \ V ub \ from the nonlcptonic 
data. The semileptonic fit gives /+(0) = 0.25 ± 0.06, so 
Eq.© yields \V ub \ nonlep = (2.9 ± 1.0) x 10~ 3 . 

Our final result for \ V ub \ is given in (|20|) . The final the- 
ory error is dominated by the lattice points, and is very 
close to their error. It will decrease with this error in 
the future. See also (^j. To go beyond the analysis here 
it will be interesting to study the additional error cor- 
relation implied by the dispersion relations when lattice 
input points are included that are closer together. 
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